A new approach to the Cramer-Rao type bound 
of the pure state model 



Keiji Matsumoto [] 
Abstract 

In this paper, new methodology - direct approach - for the deter- 
mination of the attainable CR type bound of the pure state model, 
is proposed and successfully applied to the wide variety of pure state 
models, for example, the 2-dimensional arbitrary model, the coherent 
model with arbitrary dimension. When the weight matrix is SLD 
Fisher information, the bound is determined for arbitrary pure state 
models. Manifestation of complex structure in the Cramer-Rao type 
bound is also discussed. 
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1 Introduction 

The quantum estimation theory deals with determination of the density 
operator of the given physical system from the data obtained in the experi- 
ment. For simplicity, it is assumed that a state belongs to a certain subset 
M. = {p{6)\6 G B C R m } of the space of the states, which is called model, 
and that the true value of the finite dimensional parameter 9 is left to be es- 
timated statistically. In this paper, we restrict ourselves to pure state model 
case, where M. is a subset of the space V\ of pure states in d-dimensional 
Hilbert space TC {d < oo). For example, Ai is a set of spin states with given 
wave function part and unknown spin part. 
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In the classical estimation (throughout the paper, 'classical estimation' 
means the estimation theory of probability distribution), the mean square 
error is often used as a measure of error of the estimate, and the Cramer-Rao 
inequality assures that the inverse of so-called Fisher information matrix is 
the tight lower bound of covariance matrices of locally unbiased estimator 
(Ref.|). 

Analogically, in the quantum estimation theory, in 1967, Helstrom showed 
that in the faithful state model, the covariance matrix is larger than or 
equal to the inverse of SLD Fisher information matrix, and that in the 1- 
dimensional faithful model, the bound is attainable Q J/J . 

On the other hand, in the multi-dimensional model, it is proved that 
there is no matrix which makes attainable lower bound of covariance matrix, 
because of non-commutative nature of quantum theory. Hence, the measure 
of the error of the estimate which is often used is TrGVg [M] , where Vg [M] 
denotes the covariance matrix of the locally unbiased measurement M at 6 
and G is a weight matrix, or an arbitrary given mxm positive symmetric real 
matrix. The infimum of is said to be attainable or achievable Cramer-Rao 
(CR) type bound of the model at 9 with weight matrix G, and to determine 
the attainable CR type bound long had been one of the main topics in this 
field, and is solved only for the several specific models, because TrGVg [M] 
is a functional of probability valued measure, or pair of infinite number of 
operators in the infinite dimensional Hilbert space. 

Yuen, Lax and Holevo found out the attainable CR type bound of the 
Gaussian state model, which is a faithful 2-dimensional model obtained by 
superposition of coherent states by Gaussian kernel[14]||]. Nagaoka and 
Hayashi calculated the attainable CR type bound of the faithful faithful 
spin-1/2 model [ 11 Fujiwara and Nagaoka determined the bound for 
the 1-dimensional pure state model and the 2-dimensional coherent model, 
which is the pure-state-limit of the Gaussian model. 

All of their works are based on a methodology, which we call indirect 
approach hereafter; First one somehow find an auxiliary bound which is not 
generally attainable and then proves it to be attained in the specific cases. 
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In the approach in this paper, called direct approach in contrast with in- 
direct approach, we reduce the problem to the minimization of the functional 
of the finite numbers of the finite dimensional vectors. 

The methodology is successfully applied to the general 2-dimensional 
pure state model, and coherent model with arbitrary dimension. These are 
relatively general category in comparison with the cases treated by other 
authors. Also, when the weight matrix is SLD Fisher information matrix, 
which will be defined in somewhere in the paper, the bound is calculated for 
arbitrary pure state models. 

As a by-product, we have rather paradoxical corollary, which asserts that 
even for 'non-commutative cases', simple measurement attains the lower 
bound. 

The paper is organized as follows. In section 2 and 3, basic concepts of 
the quantum estimation theory are introduced. In section 4, the commuting 
theorem, which plays key role in the foundation of the direct approach, is 
presented and is applyed to the characterization the quasi-classical model, in 
which non-commutative nature of the theory is not apparent. We formulate 
the problem in the non-quasi-classical models in section 5. Our new method- 
ology, direct approach, is introduced in section 6 and 7, and is applied to 
the 2-dimensional pure state model and the coherent model in section 8 and 
11 respectively. In section 9, we consider informational correlation between 
the parameters, and the attainable CR type bound for the direct sum of 
the models. The manifestation of the quantum structure, together with the 
minimization of the minimum of Tr J s (9)Vg[M], is discussed in section 10. 

2 Locally unbiased measurement 

Let <r(R m ) be a a- field in the space R m . Whatever measuring apparatus 
is used to produce the estimate 9 of the true value of the parameter 9, the 
probability that the estimate 9 lie in a particular measurable set B in R m 
will be given by 

Pi{9 e B\9} = trp{9)M(B) (1) 
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when 9 represents the true value of parameter. Here M is a mapping of a 
measurable set B E <r(R m ) to non-negative Hermitian operators on H, such 
that 

M(4>) = 0,M(R m ) = /, 

oo oo 

M(\jB l ) = Y / M(B l ) (B i nB j = <f>,i^j), (2) 
i=i i=i 

(see Ref.|7|,p.53 and Ref.[]|],p.50.). M is called a generalized measurement 
or measurement, because there is a corresponding measuring apparatus to 



any M satisfying (j2|) [fT^1 |lS]. A measurement E is said to be simple if E is 
projection valued. 

A generalized measurement M is called an unbiased measurement in the 
model M, if E e [M] = 9 holds for all 9 E G, i.e., 

| l^'trpW Af ((cW) = (j = 1, • • • , m). (3) 
Differentiation yields 

J§hT^-M(d9)=6i, (j,k = l,...,m). (4) 

If (|3|) and (^) hold at a some 0, M is said to be locally unbiased at 9. 
Obviously, M is unbiased iff M is locally unbiased at every 9 E 0. 

As a measure of error of a locally unbiased measurement M, we employ 
the covariance matrix with respect to M at the state p$, V$[M] = [v J e k ] E 
R mxm , where 

vl k = f(9 j -9 j )(9 k -9 k )tvp(9)M(d9). (5) 



We often abbreviated notation V[M] for Vg[M] when it is not confusing. 
The problem treated in this note is to find a lower bound for Vg[A/]. 

Only locally unbiased measurements are treated from now on, because 
of the following reason. Given N copies of the system, we apply a proper 
measurement to the the first pN copies, and the true value of parameter 
is known to lie in certain e-ball centered at 9$ with the probability ~ 1 — 
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e -a/e 2 N _ Therefore, applying the 'best' locally unbiased measurement at 0$ 
to the (1 —p)N copies, we can achieve the efficiency arbitrarily close to that 
of the 'best' locally unbiased measurement at 9, in the sense of the first 
order asymptotics. 

3 CR bound by SLD Fisher information matrix 

In 1995, Fujiwara and Nagaoka (3| defined SLD Fisher information for pure 
state models. Here, we try another definition which is adequate for our 
direct approach. 

Analogically to the classical estimation theory, in the quantum estima- 
tion theory, we have the following SLD CR inequality, which is proved for 
the exact state model by Helstrom |||7|], and is proved for the pure state 
model by Fujiwara and Nagaoka M : 

V e [M) > (J s (0))-\ (6) 

i.e., Vq[M] - (J s (9))' 1 is non-negative definite. Here J s {9), called SLD 
Fisher information matrix, is defined by 

j s (e) = [Re(k(e)\m)}, 

where the notations \h(9)) (i = 1, ...,m) are defined afterward. 

The inequality (|6|) is of special interest, because J s ~ 1 (9) is the one of 
the best bounds in the sense of the following theorem, which will be proved 
in the section ||. 

Theorem 1 Letting A be a real hermitian matrix which is larger than J 5-1 , 
that is, A > J s ~ l , there exists such an unbiased estimator M that V[M\ is 
not smaller than A. 

To define the notations \h{9)) (i = l,...,m) and to prove the SLD CR in- 
equality, we introduce some basic notations. Ti is a set of vectors with unit 
length, 

n = {\4>)\\4>)eH,{<j ) \4>) = i}. 
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V\ denotes the totality of density operators of pure states in H. A map ir 
from H to V\ is defined by 

For the manifold M = {\</>(6)) \ 9 E 6 C R m } in H, ir(N) is defined to be a 
manifold in V\ such that 

= {p{9) | p{9) = tt(|0(0)», \<P{6)) G AA}. 

Through out the paper, we only treat with the pure state model A4 which 
writes M = tt(M) for a manifold J\f in V\. 

The horizontal lift \lx) of a tangent vector X E 7^(0) (.M) to | </>(#)), is an 
element of TC which satisfies 

Xp{0) = \[\lx){<t>{0)\ + \m){lx\), (7) 

and 

{i x \m) = o- (s) 

Here, X in the left hand side (|7|) of is to be understood as a differential 
operator. We use the symbol \h(9)) to denote a horizontal lift of di E 

Notice that spann{\li) \ i = l,...,m} is a representation of T p ^(Ai) 
because of unique existence of the horizontal lift to | </>(#)) which is proved 
as follows. Application of a differential operator X to the both sides of 
p(Q) = p2(0) yie i ds 

Xp{9) = {Xp[6))W))W)\ + \m)(m\(Xp(9)) (9) 

and therefore \lx) is given by (^Xp(9))\<p(9)}. Actually, taking trace of 
both sides of @, it is shown that (X p(9))\4>{9)} satisfies (||). To prove the 
uniqueness, it suffices to show that |Z) = if (l\cp(9)) = and 

o = \i){<Ko)\ + \<Km\ (io) 
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holds true. Multiplication of \4>(9)) to the both sides of ( |l0| ) proves the 
statement. 

Fujiwara and Nagaoka defined SLD Fisher information matrix J s (9) by 
using the symmetric logarithmic derivative (SLD) of the parameter 9 l is a 
hermitian matrix Lf(9) which satisfies 

aiP(e) = ±{Lf(e)p(6) + P {o)Lf(o)). (11) 

Using SLD, the horizontal lift of $ to T\ m) (H) writes |P(0)) = Lf (0)|0(0)). 
J s (0) is called SLD Fisher information matrix because J s (0) writes 

J 5 (0) = [Retrp(0)Lf(0)Lf(0)]. (12) 

SLD defined by ( pd| ) has the arbitrariness which corresponds to the kernel of 
p(9), and Fujiwara and Nagaoka || showed that J (0) is uniquely defined 
regardless this arbitrariness. Notice that in our framework, uniqueness of 
SLD Fisher information matrix is trivial. 

We define estimation vector \ x l [M, |</>(0))]) of the parameter 0* by a mea- 
surement M at |0(0)), by 

|^[M,|^(0))]) = j^-^MWI^)}. 

An estimation vector |x*[M, |0(0))]) is said to be locally unbiased iff M is 
locally unbiased. The local unbiasedness conditions for estimating vectors 
writes 

(x l [MM0))]\m) = 0, (13) 
Re(x i [M,|</»(0))]|/ j (0)) = 6 i j (i,j = l,...,m). (14) 

Often, we omit the argument in \lj(9)), |0(0)),p(0), and J (0) and denote 
them simply by |^), \<j>),p, J s ■ Also, |a;*[M, |0(0))]) is denoted simply by \x l ) 
with the arguments M, \<j>(9)) left out, so far as no confusion is caused. 
We denote the ordered pair of vectors 
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and 

[|zi),K2>,-lC>] 

by X and L respectively. Then, the unbiasedness conditions ( |LH ) writes 

ReX*l = Re[(x i \l j )}= I m , (15) 

where I m is the mxm unit matrix. The SLD Fisher information matrix J s 
writes 

J s = ReL*L 

The imaginary part of L*L is denoted by J. 

Now, we are in the position to prove SLD CR inequality. 

Lemma 1 Following two inequalities are valid: 

V[M] > ReX*X. (16) 
V[M] > X*X. (17) 



Lemma 2 



ReX*X > J 



(18) 



holds. The equality is valid iff 



or, equivalently, 



x = u 8 - 1 = 



.., m 



(19) 
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They are proved in almost the same manner as the strictly positive case 
(see Ref.f| p. 88 and p.274 respectively). Lemmas [JD !ead to the SLD CR 
inequality @. 

Theorem 2 (Fujiwara and NagaokaJ^j) SLD Fisher information gives a 
lower bound of covariance matrix of an unbiased measurement, i.e., (||) holds 
true. 

The SLD CR inequality @ looks quite analogical to CR inequality in 
classical estimation theory. However, as is found out in the next section, the 
equality does not generally establish. 

4 The commuting theorem and the quasi-classical 
model 

In this section, the necessary and sufficient condition for the equality in the 
SLD CR inequality to establish is studied. Fujiwara has proved the following 
theorem 

Theorem 3 (Fujiwara J^j) The equality in the SLD CR inequality estab- 
lishes iff SLDs {Lf \i = 1, ...,m} can be chosen so that 

[Lf,Lf] = 0, = l...,m). 



We prove another necessary and sufficient condition which is much easier to 
check for given models, by use of the following commuting theorem, which 
plays key role in our direct approach. 

Theorem 4 If there exists a unbiased measurement M such that 



\x>) = j{F-P)M{d0)\<t>), 

V[M] = ReX*X, (20) 
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then, 



ImX*X = (21) 

holds true. On the other hand, if ( pC| ) holds true, then there exists such a 
simple, or projection valued, unbiased measurement E that (^) holds and 

E{{6 K })E({6 K }) = E({9 K }), 
E({§ }) = Eq, 

e(b™/[J{O k })=0, (22) 

for some {9 K \9 K G H m ,K = 0, ...,m + 1}, where Eq is a projection onto 
orthogonal complement subspace of spanc{X-} . 

Proof If (|2^) holds, inequality ((I?]) in lemma [l| leads to 

ReX*X > X*X, 

or 

> fImX*X, 

which implies ImX*X = 0. 

Conversely, Let us assume that (^T|) holds true. Applying Schmidt's 
orthogonalization to {\(f>), {x 1 ), |x m )} and normalizing the product of or- 
thogonalization, we obtain the orthonormal system {\b l ) \ i = 1, ...,m + 1} 
of vectors such that, 

m+l 

k l )= E A il^)' 3A|eR,i = l,.., m ,i = l,.., m + l. 
i=i 

Letting O = [o*-] be a (m + 1) x (m + 1) real orthogonal matrix such that 

m+l 
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and denoting Y^f=i by W l )i we have 



x 



m+l m+l 

E A ;E ^ ,fc > 
j=i fc=i 

m+l /m+l \ 

E E l**> 

fc=i \i=i / 



m+l v^m+1 \ i fc 

Therefore, noticing that the system {\b H ) \ i = l,...,m + 1} of vectors is 
orthonormal, we obtain an unbiased measurement which satisfies (22) as 
follows: 

e *= WW) ' « = L + 

O = O, 

£(R}) = |6" s )(6" s |, « = l,...,m + l, 

m+l 

E(0 o )=I H - \b' K )(b' K \, 

K=l 

E(R m /{6 ,..J m+ i}) = 0. 
Here, /-^ is the identity in H. □ 



Theorem 5 The equality in the SLD CR inequality establishes iff 

ImL*L = (23) 

{lj\h) is real for any When the equality establishes, that bound is 

achieved by a simple measurement, i.e., a projection valued measurement. 

Proof If the equality establishes, by virtue of lemma [lf]2|, we have (|2l|) 
and (|19|), which lead directly to (|23]). 

Conversely, if Jm(lj\lk) = for any j, k, by virtue of commuting theorem, 
there exists such a simple measurement E that 

YH jS ~ 1 ) j '%) = [0-e)E(d§)\<t>). 
k J 



li 



Elementary calculations show that the covariance matrix of this measure- 
ment equals J 5-1 . □ Our theorem is equivalent to Fujiwara's 

one, because by virtue of commuting theorem, (lj\h) is real iff there exist 
such SLDs that Lf and Lj commute for any i,j. However, our condition is 
much easier to be checked, because to check Fujiwara's condition, you must 
calculate all the possible SLDs, for the SLD is not unique. In addition, SLD 
is much harder to calculate than horizontal lift. When the model has only 
one -dimensional, we have the following corollary of theorem [5] 

Corollary 6 when a manifold Ai is one-dimensional, the inverse of SLD 
Fisher information matrix is always attainable by a simple measurement. 

Remark Often, a model is defined by an initial state and generators, 

N = {\4>(0)) | \(f>(0o)) = |^o); di\4>(6)) = iHi(6)\(j){6)), 6> e 6 C R" 1 }, 
M = tt(M). 

Then, (lj\k) is real iff (c/)(6)\lHi(0), Hj(0)}\(j)(0)) = 0, which is equivalent to 
the existence of generators which commute with each other, [Hi(9), Hj(9)] = 
by virtue of the commuting theorem. 

Putting the remark and the Fujiwara's theorem together, we may metaphor- 
ically say that the equality in the inverse of SLD Fisher information matrix 
is attainable iff any two parameters 'commute' at 9. Throughout the pa- 
per, we say that a manifold Ai is quasi- classical at 9 iff {lj\h) is real at 9. 
The following remark describes another 'classical' aspect of the condition 
Imiljlk) = 0. 

Example When the model Ai is given by 

M = {p{0) | p(9) = ir(\4>(9))), \4>(9))is an element of real Hilbert space}, 

the model is quasi-classical at any point in A4. 

As is illustrated in this example, when th model Ai is quasi-classical at 
9q, a state vector \4>(9)) behaves like an element of real Hilbert space around 
#o, and the state vector's phase parts don't change around #o at all. 
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5 Non-quasi-classical cases 



As was concluded, the equality in the SLD CR inequality establishes only 
when the model is quasi-classical, and there is not any better bound than 
the inverse of SLD Fisher information matrix, as in theorem |l|, which is 
straightforwardly derived from the following lemma, which is proved in the 
appendix A. 

Lemma 3 For any i, 

inf {[V&[M]]m I M is locally unbiased at 6} = \j s ' 1 (9) 



In general case, therefore, we must give up to find a matrix which makes 
attainable lower bound of V[M], and instead, we try to determine 



CR(9,G,M) = mi{TrGV | V € V d {M)} 



(24) 



for an arbitrary nonnegative symmetric real matrix G, where Vo{M) (, or in 
short, V,) is the region of the map Vg[*] from unbiased estimators to m x m 
real positive symmetric matrices. CR(#,G, Ai) is the attainable CR type 
bound, and we often use abbreviated notations such as CR(#, G), CR(G) . 

To make the estimational meaning of ( |24] ) clear, let us restrict ourselves 
to the case when G is diag(gi , g2---,9m)- Then, the attainable CR type bound 
is nothing but the weighed sum of the covariance of the estimation of 6 l . If 
one needs to know, for example, 9 1 more precisely than other parameters, 
then he set g\ larger than any other g^, and choose a measurement which 
achieves the attainable CR type bound. 

Notice that 



> 



inf | X>[ty[M]]i 

m 

£inf{^[^[Ml]« 

i=\ 

Eft V S ~ X H^ 



M is locally unbiased at 6 j 
M is locally unbiased at 9} 
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holds true by virtue of the lemma [31 and that the equality in the first inequal- 
ity does not always establish, implying that in the simultaneous estimation of 
different parameters, there is information losses because of non-commutative 
nature of the quantum mechanics. 

Another proper alternative of the classical Fisher information matrix is 
a set inf Vq(M) of symmetric real matrices, where the notation inf is defined 
as follows. Let us define 

IbX = {A | A is real and symmetric, V-B < A, B G X}, 
ubX = {A | A is real and symmetric, \/B > A, B E X}, 

where X is a set of real symmetric matrices, and we define inf V by 

infV = lbVnub(lbV). 

Then, we have the following lemma. 

Lemma 4 infV is a subset of the boundary bdV ofV. 

This lemma is a straightforward consequence of the following lemma, which 
is proved in the appendix B. 

Lemma 5 If V is an element of V, then V + Vq is also an element of V, 
where Vq is an arbitrary real nonnegative symmetric matrix. 

Because of lemma |], it is of interest to determine the boundary bdV. bdV 
is turned out to be a subset of V such that V = CR(G) for a weight matrix 
G, because of lemma ||, and lemmas ^-0. 

Lemma 6 V is convex. 

Proof Let M and iV be an unbiased estimator. Because 

XV[M] + (1 - X)V[N] = V[\M + (1 - X)N] 

holds true and AM + (1 — A)A^ is an unbiased estimator, we have the lemma. 
□ 
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Lemma 7 V is closed. 



Lemma ^ will be proved in the appendix C. 



If a model A4 has smaller value of the attainable CR type bound at 9 
than another model M at 9' has, the Vg(A4) of is located in the 'lower part' 
of Sym{m) compared with that of Vqi(J\T). 

6 The reduction theorem and the direct approach 



Any generalized measurement M in TC can be dilated to a simple measure- 
ment E in a larger Hilbert space TC' D TC, so that 



will hold, where P is the projection from TC' onto TC. 

Naimark's theorem, mixed with commuting theorem, leads to the following 
reduction theorem, which is essential to our direct approach. 

Theorem 8 Let M be a m-dimensional manifold in V\, and Bq be a system 
| i = 1, ...,m} of vectors in 2m + 1-dimensional Hilbert space 7i' e 
such that 



for any i,j. Then, for any locally unbiased estimator M at 9 in 7i, there is 
a simple measurement E in Ti' e such that 'locally unbiasedness' is satisfied, 



Theorem 7 (Naimark's theorem, see Ref. Q, pp. 64-68.) 



M (B) = PE(B)P 



(25) 



mi'j) 




(26) 



(x l W) = o, 
Re{x%) = 8) 




i,j = 1, ... 



m) 



(27) 
(28) 
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and that the 'covariance matrix' V[E] of E equals V[M], 



V[M] = V[E] 



{e 1 - ei)txpE(de) 



(29) 



Proof For any locally unbiased measurement M, there exists a Hilbert 
space TLm and a simple measurement Em in T~Lm which satisfies ( f25|) by 
virtue of Naimark's theorem. Note that Em is also locally unbiased. Let 
\y l ) £ TCm denote the estimation vector of Q l by Em, that is, 

\ y i)= j{fr-&)E M (dB)\<f>'). 

Mapping spanc{\4>), |Zj), \y l ) \ i = l,...,m} isometrically onto so that 
{\<p),\li) | z = l...m} are mapped to {|</>'),|^) | % = l,...,m} , we denote the 
images of {\y l ) \ i = 1, m} by {\x l ) \ i = 1, m}. 

Then, by virtue of the commuting theorem, we can construct a simple 
measurement E in 7i' g satisfying the equations fl2"6|) - ( 2£ ) . □ 



The reduction theorem shows that V is identical with the set of matrices 

V = ReX*X 

such that 

\x i )eH' e \{W)}{i = l,...,m), 

where 7i' e \{\4>')} denotes the orthogonal complement subspace of Tt' e , and 
that ( |2~l| ) and ( |2"8| ) are satisfied. Now, the problem is simplified to the large 
extent, because we only need to treat with vectors {\x l ) \ i = 1, ...,m} in 
finite dimensional Hilbert space H.' e instead of measurements, or operator 
valued measures. 

We conclude this section with a corollary of reduction theorem, which is 
rather counter-intuitive because historically, non-projection- valued measure- 
ment is introduced to describe simultaneous measurements of non-commuting 
observables. 
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Corollary 9 When the dimension of TC is larger than or equal to 2m + 1, 
for any unbiased measurement M in TC, there is a simple measurement E 
in TC which has the same covariance matrix as that of M . 

Proof Chose {\l[) \ i = l,...,m} to be \ i = l,...,m}. □ Especially, if 

TC is infinite dimensional, as is the space of wave functions, the assumption 
of the corollary is always satisfied. 

7 Lagrange's method of indeterminate coefficients 
in the pure state estimation theory 

Now, we apply our direct approach to the problem presented in the section 
||, or the minimization of the functional TrGReX*X of vectors in TC'g. One 
of most straightforward approaches to this problem is Langrange's indeter- 
minate coefficients method. First, denoting an ordered pair {\l[ \ i = 1,...,} 
of vectors in TC' e also by L, the symbol which is used also for an ordered pair 
{\li | i = 1, } of vectors in TC, we define a function Lag(X) by 

Lag{X) = ReTrX*XG - 2Tr((ReX*L - J m )2) - TrImX*XA, (30) 

where 3, A are real matrices whose components are Langrange's indetermi- 
nate coefficients. Here, A can be chosen to be antisymmetric, for 

TrImX*XA = TrImX*X(A - A T )/2 

holds true and only antisymmetric part of A appears in (|3C|). 

From here, we follow the routine of Langrange's method of indeterminate 
coefficients. Differentiating L(X + e<5X) with respect to e and substituting 
into e in the derivative, we get 

ReTr(£X*(2XG - 2LH - 2iXA)) = 0. 

Because 5X is arbitrary, 

X(G - iA) = LH (31) 
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is induced. 

Multipling X* to both sides of (|3l|), the real part of the outcomming 
equation, together with (|l5|), yields 

S = ReX*XG = VG. (32) 

into ([H]) , we obtain 

X(G - iA) = LFG. (33) 

In this paper, we solve fli"5|), (|2l|), (^) and V = ReX*X with respect to 
X, real symmetric matrix V and real antisymmetric matrix A, for the variety 
of pure state models. However, the general solution is still far out of our 
reach. 



Substituting (32 



8 The model with two parameters 

In this section, we determine the boundary of the set V in the case of the 
2-dimensional model. 

The equation(|33|), mixed with (|2l|), leads to 



(G-iA)V(G-iA) = GVl*lVG. (34) 
whose real part and imaginary part are 

GVG - AVA = GV J S VG, (35) 

and 

GVA + AVG = -GVJVG, (36) 

where J denotes ImL*L, respectively. 

As is proved in the following, when the matrix G is strictly positive, 



( p4|) is equivalent to the existence of X which satisfies ( 15 ) , (|2lD , ((JJ) , and 



V = ReX*X. If V and A satisfying (M) exist, X which satisfies (S3) and (EB) 
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is given by X = UV 1 / 2 , where U is such a 2m + 1 X to complex matrix that 
U*U = I m . X = UV 1 / 2 also satisfies (|l5|), because 

VG = ReX*lVG 

is obtained by multipling X* to and taking real part of the both sides of 

(ED- 

Hence, if G is strictly positive, our task is to solve ( |35"1 ) and (^) for real 
positive symmetric matrix V and real antisymmetric matrix A. When G 



is not strictly positive, after solving (35) and (|36|), we must check whether 
there exists such X which satisfies (|l5"|), (21) and V = ReX*X. 

Throughout this section, we parameterize the model so that J s is equal 
to the identity matrix I m . Given an arbitrary coordinate system {6 l \ i = 



l,...,m}, such a coordinate system {&'' 
following coordinate transform: 



l,...,m} is obtained by the 



e' l = Y,l(J s ) 1/2 k° J d = h-,m). 

By this coordinate transform, V is transformed as: 

V' = {J s f/ 2 V(J s f/ 2 . 



(37) 



(38) 



If the result in the originally given coordinate is needed, one only needs to 



transform the result in the coordinate system {9' 



1, to} using 



in the converse way. 

So far, we have not assumed dimTW = 2. When dim.M = 2, covariance 
matrices are included in the space Sym{2) of 2 x 2 symmetric matrices which 
is parameterized by x, y, and z, where 



V 



V 



z + x y 
y z — x 



Sym{2) 

Before tackling the equations (|35|) and (ptf), three useful facts about this 



parameterization are noted. First, letting A is a symmetric real matrix 
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which is represented by (A x ,A y ,A z ) in the (x, y, z)-space, the set C + (A) of 
all matrices larger than A is 

C+{A) = {(x, y, z)\(z- A z f - (x - A x ) 2 - (y - A y f > 0, z > A z }, 

that is, inside of a upside-down corn with its vertex at A = (A x , A y , A z ). 
Hence, V is a subset of C+(J m ), or inside of a upside-down corn with its 
vertex at (0, 0, 1) because of the SLD CR inequality. When the model M. is 
classical at 9, V coincides with C + (I m ). 

Second, an action of rotation matrix Rq to V such that RgVRj , where 

cos 6 — sin 6 
Re = i 
sin 6 cos 6 

corresponds to the rotation in the (x, y, z)-space around z-axis by the angle 
29. 

Third, we have the following lemma. 

Lemma 8 V is rotationally symmetric around z-axis, if M is parameterized 
so that J s writes the unit matrix I m . 

Proof The necessary and sufficient condition for V to have rotational 
symmetry around z-axis is the existence of a 2m + 1 by m complex matrix 
Y satisfying @, @ and 

yd R e (ReX*X) Rj = R e X*XR% = Y*Y, (39) 

for any given ordered pair X of vectors which satisfies (|l5|) and 

On the other hand, because of L*L = I m + iJ, elementary calculation 
shows 

L*L = R e l*lRj, 
or equivalently, for some unitary transform in 7i' e \ {| </>)}, 

l*U = R e l*, 
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which leads, together with (15), to 

Rel*UX 

Therefore, 



satisfies 



Y = UXRj 
and we have the lemma. 



□ 



Now, the boundary of the intersection V of V and zx-plain is to be 
calculated, because V is obtained by rotating V around z-axis, by virtue of 
this lemma. bdV is obtained as the totality of the matrix V = ReX*X which 
satisfies (|l5|) , (|2l|), and (|33|), for a diagonal real nonnegative matrix G. 

Let us begin with the case where a diagonal matrix G is positive definite. 
In this case, we only need to deal with (5^) and (|36|). Let 



J 



-p 

(3 



and 



A 



" —A 


,G = 


' 1 


" 




' u " 








,v = 




. A 







9 _ 




v 



where g, v, and u are positive real real numbers. Note that 

\fi\<l 

holds, because L*L = I m + J is nonnegative definite. Then, ( |35"1) and 
writes 



u + v\ 2 - u= 0, 
vg 2 + u\ 2 - v 2 g 2 = 0, 

vgX + uX + uvfig = 0. (40) 

The necessary and sufficient condition for A and positive g to exist is, after 
some calculations, 



y/u-l + y/v-1- \P\y/uv = 0. 



(41) 
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Note that u and v are larger than or equal to 1, because V > J s 1 = I m . 
Substitution of u = z + x and v = z — x into (P) and some calculation leads 
to 

\/3\ ,J(z + x - l)(z - x - 1) ± yjl - f3 2 (y~z + x - 1 + Vz - x - l) = \p\. 

It is easily shown that the lower sign in the equation corresponds to the set 
of stationary points, and 



(z + x- l)(z -x-l)±Jl-(3 2 (Vz + x-l + Vz-x-l) = \p\ 



(42) 

gives a part of bdV. In fl42]), x takes value ranging from — /3 2 /(l — 1 ) to 
f3 2 /(l — f3 2 ) if is smaller than 1. When \ j3\ = 1, x varies from — oo to oo. 
This restriction on the range of x comes from the positivity of z — x — 1 and 
z + x — 1. 
When 



G 



' 1 


" 


, or G = 


' 


" 










_ 


_ 







1 



(43) 



we must treat the case of = 1 and the case of < 1 differently. In the 
case of |/3 1 = 1, there exists no 2m x m complex matrix X which satisfies 
V = ReX*X, (0) and @). On the other hand, if \0\ < 1, such complex 
matrix X always exists and V = ReX*X is given by, in terms of (x, y, z), 



—x + 1, x < 
or 

x + 1, x > 



P 2 



l-P 2 



ft 2 



l-(3 2 



(44) 



Because any element on the line (fl4|), if x 7^ ±/3 2 / (1— /3 2 ), has an element 
of V which is smaller than itsself, (|42] ) the intersection of inf V and zx-plane, 
where 



(3 2 /(l-f3 2 )<x<(3 2 /(l-(3 2 ). 
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Figure 1: Two stationary lines; 

(i) \P\y/{z + X-1){Z-X-1) + ^Jl^W (y/z + X-l + y/z-X-^ 

(ii) \0\y/(z + x-l){z-x-l)- y/T^W (^z + x-l + ^z-x-l) 
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The intersection of z-axis and bdV gives 

CR(J S ) = . It> , (45) 

where the equality holds in any parameterization of the model A4 . 

In arbitrary parameterization of the model M, with help of ( |42"D and 
(|38p, inf(V) is obtained as, for \/3\ > 0, 

det yfi(yj + yfojfP) - 1 Tr^(V) = 1 (\(3\ > 0), 



TryJV(V) = {p = 0), (46) 

where 

V(V) = VJsv^/Js - I m . 
Slight look at the equations (pq ) leads to the following theorem. 
Theorem 10 In the 2- dimensional model, if 

\p(e,M)\>\p(e',M% 

J S (6,M) = J S (6',M'), 
then the Vg(M) is a subset ofVg>(M'). 

The equations ( f46| ) and tedious but elementary calculations shows the fol- 
lowing theorem. 

Theorem 11 In the 2- dimensional model, if 

\P(9,M)\ = \(3(6',M')\, 
J S (6,M) < J S (8',M'), 

then the Ve(M) is a subset ofVg'(M'). 

By virtue of these theorems, \/3\ can be seen as a measure of 'uncertainty' 
between the two parameters. Two extreme cases are worthy of special at- 
tention; When \/3\ = 0, the model M is classical at 6 and V is maximum. 
On the other hand, if |/3| = 1 , V is minimum and uncertainty between 9 1 
and 9 2 is maximum. In the latter case, we say that the model is coherent at 

e. 
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Example (spin rotation model) We define spin rotation model Q by 
M,, m = 7r(A0> 

M s ,m = {\4>{0)) | \4>(0)) = T(0)\s,m), < 6 1 < ir, < 2 < 2vr} , 

T(0) = exp (# (sin 6 2 S X - cos fl 2 ^) , (47) 
where S x , S y , S z are spin operators, and \s,m) is defined by, 
Sz\j, m ) = hrn\s,m), 

(Sl + S 2 + S 2 ) \s,m) = h 2 s{s + l)\s,m). 

s takes value of half integers, and m is a half integer such that —j<m<j. 
Then after tedious calculations, we obtain 

rn 

0(0,M .,m) = -2T 2- 

s z + s — m z 

If m = as, where a < 1 is a constant, f3(0,M s ,m) tends to zero as s — > oo, 
and the model M. s ,m becomes quasi-classical. However, if m = s, the model 
M. s ,m is coherent for any s. 

Example (shifted number state model) shifted number state model, which 
has four parameters, is defined by 

M n = ir{Nn), 

N n = [\<t>{d)) | \<f>(6)) = D{6)\n), G R 2 } , 

where letting P, X be the momentum operator and the position operator 
respectively, 

D(8) = exp 1 (-6 l X + 6 2 P) , 
and \n) is the nth eigenstate of the harmonic oscillator, 

H = --P 2 + -X 2 . 
2 2 

After some calculations, we have 

-1 



0(0, M n ) 



2n+l 

As n tends to infinity, (3(9, Ai n ) goes to and the model becomes quasi- 
classical. 
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9 Informational exclusiveness and independence, 
and direct sum of the models 



In a m-dimensional model M, we say parameter 9 l and 3 are information- 
ally independent at 6>o, iff 

Re(li\lj)\o = o = Im(Zj|Zj)|0 = o = 0, 

because, if the equation holds true, letting the submodels M(1\9q), •M.(2\6q) 
and M(l,2\9 ) of M, be 

M(i\e ) = {p(e)\e = (e\el...,e^),6 1 eK}, 
M(2\e ) = {p(9)\e = (e 1 Q ,e 2 ,...,e^),e 2 eK}, 

M(l,2\9 ) = {p(6)\8 = (e\9 2 ,9l...,9Z l ),(6\9 2 )£B 2 }, (48) 
the following equality establishes: 

CR(9 , diag( gi , g 2 ), M(l, 2\6 )) = CR{9 , diag( gi ), M(l\9 )) + CR(9 , diag(g 2 ), M(2\9 )), 

which means that in the simultaneous estimation of the parameter (9 1 ,9 2 ), 
both of the parameters can be estimated without the loss of information 
compared with the estimation of each parameters. 
On the other hand, iff 

Re{h\l 2 )\e=e = 0, (49) 
and A^(l,2|6*o) is coherent, or equivalently, 

lm(h\l 2 )\ e =e = mh)(h\h)) 1/2 \ e=eo 

hold true, we say the parameters are informationally exclusive at 9o, because 
of the following theorem. 

Theorem 12 Let 9 1 and 9 2 be informationally exclusive parameters at 9o, 
and M' a measurement which takes value in R 2 and satisfies local unbiased- 
ness condition about 9 1 at 9q, 

J 9 1 TT P {9 )M'{d8) = Ol 



! 



M'{d9) = (50) 

0=00 
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If the measurement M' satisfies i.e., 

(0 1 -9 1 ) 2 Ttp(9 )M'(d9) = CR(diag(l,0),9 ) 



S-l 



J 



M' can extract no information about 9 2 from the system, i.e., 



11 



VBcR 2 TV M\B) 



dp 



o, 



=e 0/ 



(51) 



(52) 



and vice versa. 



Proof Let £ be a Naimark's dilation of the measurement M' , and de- 
compose the estimation vector \x) = \x[E, |</>($o))]) of E as 

\x)=z\<f>(6 o ))+w\h(0 Q )) + \^), 

where \ip) is orthogonal to both of \4>(9o)) and \li(6o))- Then, local unbi- 

i li 



asedness condition leads to z = and w 



J 



s-i 



must be the 



zero vector for M' to achieve the the equality (51), because the variance of 
M' writes 



(9 1 - 9l) 2 Tr p(9 )M'(dB) = (x\x) 



J 



s-i 



ii 



Using the fact that by virtue of informational exclusiveness, 1/2(^0)) 
writes 

IM0o)>=io|Zi(0 o )>, 

where i is the imaginary unit and a a real number, we can check the equality 
(|52"D by the following calculations: 



Tr M ' (B) w 



=0o. 



Rb{</>(0 o )\E(B)\1 2 ) 
-alm{<t>{0 o )\E(B)\l x ) 



J" 



11 



Im {<P(9 )\E(B) / (6 L - e L )E{dO)\<t>{0 )) 



2N 



□ 



Fujiwara and Nagaoka Q showed that in the 2-dimensional model with 
the informationally exclusive parameters, the best strategy for the estima- 
tion is alternative application of the best measurement for each parameter 



to the system. This fact is quite natural in the light of theorem 12 
For the submodels 

Mi = M(l,2, ...,rai | 6> ), M 2 = M(mi,mi + 1, ...,m\ 9 ) 



of M, which are defined almost in the same way as the definition (48) of 
M(l, 2 | #o); we sa Y that M is the sum of M and M at 6q, and express the 
fact by the notation, 

M\o =M 1 ®M 2 \e - 
Throughout the section, m — mi is denoted by m 2 . 

Lemma 9 // any "parameter of Mi is informationally independent of any 
parameter of M 2 at 9q, and the weight matrix G writes 



G 



Gi 
G 2 



then 



CR(G,9 ,M) = CR(Gi,e ,Mi) + CR(G 2 ,e ,M2) 



When the premise of the lemma is satisfied, Mi and M 2 are said to be 
informationally independent at #o- 

Proof Let M be a locally unbiased measurement in M, and define the 
measurements Mj (j = 1, 2) in R m J (j = 1, 2) by 

Mi{B) = M(B x R™ 2 ) (B G R mi ), 
M 2 (B') = M(R mi x B') (B 1 G R m2 ), 
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respectively. Then, the measurement Mj (J = 1,2) is locally unbiased in 
Aij (j = 1,2), respectively. 
Therefore, we have 

inf {TrGV[M] | M is locally unbiased in M] 
= inf {TrG\V[Mi] \ M is locally unbiased in M) 

+ inf {TrG 2 V[M 2 ] \ M is locally unbiased in M} 
> inf { TvGiV[M'] \ M' is locally unbiased in Mi} 

+ M{TrG 2 V[M"] | M" is locally unbiased in M 2 ) 



or its equivalence, 

CR(G,M) > CR(Gi,7Wi) + CR(G 2 ,M 2 ). 

Because .Mi and M 2 are informationally independent, L for A4 writes 

Li 
L 2 



(53) 



in the appropriate coordinate, where l_i = |^ mi )], and L 2 

In that coordinate, let us write X as 

Xl X12 



l mi+l/i 



Then 



ReX*U 



X 



I mi , lmX*Xi 



X21 X2 



0(* = 1,2), Xi 2 =X 2 i = 0, 



is a sufficient condition for the measurements corresponding to X to be 
locally unbiased. Therefore, we have 

CR{G,M) 

= inf { TrGX*X | ReX*L = J m , ImX*X = } 



< inf j^TTrG^Xi 



ReX*L| 



0, (i = l,2) 



= CR(G 1 ,M 1 ) + CR{G 2 ,M 2 ), 
which, mixed with (|5^) leads to the lemma. 



□ 
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10 Manifestation of complex structure 

It is worthy of notice that \/3\, which was shown to be a good index of 
'uncertainty' in the case of the 2-dimensional model, is deeply related to the 
natural complex structure in V\. 

Let us define the linear transform D in T p (M) as follows; First, multiply 
the imaginary unit i to \lx)- I n general, however, i\lx) is not an element of 
span^L, and does not represent any of vectors in T p (Ai). Hence, we project 
i\lx) onto span-Ri. with respect to the inner product Re(*|*), and the image 
by 7T* of the product of the projection is defined to be Dl G T p {M) , where 
7T* is the differential map of it. 

By elementary linear algebra, it is shown that the matrix which corre- 
sponds to D is J S ~ 1 J, and that, in the 2-dimensional model, its eigenvalues 
are ±if5. 



i\l) G spann{L, iL} 
I project 
| to) G span-RY. 

I 7T* 

DI G T p (M) 

The definition of the map D naturally leads to the following theorems. 

Theorem 13 The absolute value of the eigenvalue o/D ; or equivalently, of 
J S ~ 1 J, is smaller than or equal to 1. 

Is the eigenvalues of the linear map D a good measure of 'uncertainty' 
in the arbitrary dimensional model? If all of eigenvalues of D vanish, as 
is shown in the section the model is quasi-classical, and 'uncertainty' 
among parameters vanishes. When eigenvalues of D do not vanish, we have 
the following theorem. 



multiplication of i 
\l) G span^L ► 

T 

horizontal lift 

T 

X G T p (M) ► 

D 
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Theorem 14 For any pure state model, 

inf{TrJ 5_1 F | V G V} 



= Tr 



2 



E 



/3 : eigenvalues of D 

The estimation theoretical meaning of minjTrJ 5 ^ 1 !/ | V G V} is hard to 
verify. However, this value remains invariant under any transform of the 
coordinate in the model Ai, and can be an index of distance between V and 
J 3 - 1 . 

Proof Because min{Tr J S ~ 1 V\ V € V} is invariant by any affine coordinate 
transform in the model M, we choose a coordinate in which J s writes I m 
and J writes 



J = 





Pi 









-Pi 



































A 










-A 


































Then, The model M is decomposed into the direct sum of the submodels 
one or two dimensional M K , 

M = ($M K , 

K 

where any two submodels M. K and M. K i are informationally independent, 
and J of a two dimensional submodel M K is 

-p K ~ 
Ac o 
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Therefore, by virtue of lemma || and the equation (fS|), we have the theorem. 
□ 

11 The coherent model 

As for the model with arbitrary dimensions, the model is said to be coherent 
at 9 iff all of the eigenvalues of (J S )~ 1 J are dbi. When the model is 2- 
dimensional, this definition of coherency reduces to |/3| = 1. The dimension 
of the coherent model is even, for the eigenvalues of J s ~ l J are of the form 
or 0. 

In this section, we determine the attainable CR type bound of the co- 
herent model. The coherent model is worthy of attention firstly because the 
coherent model is 'the maximal uncertainty' model, secondly because there 
are several physically important models which are coherent. 

The definition of the map D leads to the following theorem. 

Lemma 10 The model M is coherent at 9 iff spann{iL} is identical to 
span-^L, or equivalently, iff spann{L,iL} is identical to span-^L. 

This lemma leads to the following lemma. 

Lemma 11 The model Ai is coherent iff the dimension of spanc\- is m/2. 

Proof First, we assume that 

dime spanc\- = m/2. (54) 

Because span^L is a m-dimensional subspace of span-R,{L, iL} whose dimen- 
sion is smaller than or equal to m because of (j5l|), we have span^{L,iL} = 
spanuL, or coherency of the model. 

Conversely, let us assume that the model is coherent. Taking an or- 
thonormal basis {ej | j = 1, ...m} of span^L such that ej +m /2 = De^, hori- 
zontal lifts {|j)| j = 1) ■■■m} of {ej\ j = 1, ...m} satisfy \j + m/2) = i\j), and 
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any element \u) of span-RL = spann{L,iL} writes 

m 

3=1 
m/2 

= ^2( a j +i a >j+(m/2))\j), 
3=1 

implying that the dimension of spancL is m/2. □ 

Fujiwara and Nagaoka f| determined the attainable CR type bound of 
the two parameter coherent model. In the following, more generally, we work 
on the bound of the coherent model with arbitrary dimension. Throughout 
the section, the weight matrix G is assumed to be strictly positive. 

When the model is coherent, ReL*X = I m or equivalently Rel_*(X — 
LJ 5_1 ) = 0, implies, by virtue of span^L = spann{i\-, L}, 

L*(X- LJ 5 " 1 ) = 0, 

or equivalently, 

L*X = L-LJ 5 - 1 = I m + UJ^ 1 . (55) 



Multiplication of L* to the both sides of fl33|), together with the equation 
(H, yields 



(J m + iJJ s - l ){G - iA) = (J s + iJ)VG. (56) 

By virtue of the coherency, both of the real part and the imaginary part of 
( |56| ) give the same equation, 

G + J J 8 - 1 A = J S VG, 

or 

VGVVG - VGJ^VG = (VGJ^JJ^VG) (g- 1 ' 2 AG- 1 ' 2 ) . (57) 
The antisymmetric part of the both hands of the equation yields 



Vgj 3 - 1 ,!,; 3 - 1 ^, g-^ag- 1 " 



o. 
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Therefore, letting «j and bi denote the eigenvalues of VGJ S 1 JJ s 1 VG 
and G _1 / 2 AG -1 / 2 respectively, we have 

Tr (Vgj^jj^Vg) (g- 1 / 2 kg- 1 / 2 ) 

= ^ciibi = ^2\ai\\bi\, 

i i 

where the last equality is valid because the left hand side of the equation 
( |57| ) is positive symmetricity virtue of the SLD CR inequality. 
On the other hand, from (|33|) or its equivalence, 

XVG (l m - iG- 1/2 AG~ 1/2 ) = LVVG, (58) 

we can deduce = 1 (i = 1, ...,m) as in the follows. 

The rank of the right hand side of (|58|) is equal to m/2 because G is 
strictly positive and 

rankL = dime span^Y. = m/2. 

On the other hand, 

rankX = dime spancX = diniR spanuX = m, 

where the second equality comes from ImX*X = and the last equality comes 
from ReX*L = I m . Therefore, the rank of the matrix I m — iG~ 1 / 2 AG~ 1 / 2 
must be m/2, and the eigenvalues of G~ l l 2 hG~ l l 2 are ±i. 
After all, we have 

min TrGV = TrGJ s_1 + TrabsGJ 5 " 1 JJ S ' 1 , 
VeV(M) 

where TrabsA means the sum of the absolute values of the eigenvalues of 
the matrix A. When the minimum is attained, the covariance matrix V is 
given by 

V = J 3 ' 1 + G- l ' 2 \G 1 / 2 J 3 ' 1 JJ s - l G l ' 2 \G- l l 2 , 

where \A\ = {AA*) 1 ' 2 . 

To check the coherency of the model, the following theorem, which is 
deduced from theorem [l^, is useful. 
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Theorem 15 the model is coherent at 8 iff 

\detJ s \ = |detJ|. 



Example (squeezed state model) Squeezed state model, which has four 
parameters, is defined by 

M = {p(z,0 | p(z,0 = \z,0(z,Z\, C}, 

where 

= D(z)S(O\0), 
D(z) = exp(za^ — ~za), 

5(0 = expQ(£ a t 2 -£a 2 ; 

Letting z = (6 1 + i8 2 )/2 1 / 2 , Q = {a + a))/2 l l 2 , and £ = 8 z e~ 2idA (0 < 9 3 , < 
9 A < 2tt), we have 



j 



cosh 26» 3 -sinh2/9 3 cos 28 A sinh 28 3 sin 2fl 4 

sinh 26> 3 sin 26» 4 cosh 26> 3 + sinh 26» 3 cos 28 4 






10 
-10 

-sinh26> 3 
sinh26> 3 








1 
sinh 2 26> 3 



Coherency of this model is easily checked by theorem 15 



|detJ s | = |detJ| = -sinh 2 26> 3 . 

Example (spin coherent model) As is pointed out by Fujiwara §, spin 
coherent modelM. s s , which is a special case of spin rotation model (|47|), is 
coherent. 
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Example (total space model) The total space model is the space of all the 
pure state V in finite dimensional Hilbert space Ti. By virtue of theorem [k], 
the coherency of the model is proved by checking that span^L is invariant 
by the multiplication of the imaginary unit i. Let \l) be a horizontal lift of a 
tangent vector at \<p). Then, i\l) is also a horizontal lift of another tangent 
vector at \<p}, because \4>) + i\l)dt is an element of 7i with unit length. 
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A proof of lemma 3 

Proof Let be a projection valued measure such that, 
„ m 
J^xE {i \dx) = [j S ~\0)] L j(&), 



and M p be an unbiased measurement at 9 such that 



x ... x 



+ ^ 0t + x + Ax 
Pi Pi 



X ... X 



{9 m }) = Pi E®{[x,x + A a 



M p (B X x ... x (R/ {0*}) x ... x B m ) = 



where p = [pj] (i = 1, ...m) is a real vector such that J2iPi = 1 an d Pi ^ 0> 
and Bi (i = 1, ...,m) are arbitrary measurable subset of R. Then, we have 
for any p, 

1 



[V e [M } 



Pi 



J s -He) 



which leads to 



inf { [Ve[M]]ii | M is locally unbiased at 9} 
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^2 Pi = i, Pi > o j 



< inf|[F [M p ]] 

On the other hand, SLD CR inequality leads to 

inf { [Ve[M]]ii | M is locally unbiased at 9} > [j 5-1 ^) 
and we have the lemma. 



□ 



B proof of lemma 5 

Proof Let M be a locally unbiased measurement, and v a = (vf, 
denote \/Voa, where a is a vector whose components are 1 or —1. Then, the 
measurement M' , which is defined by 

M' f jj([oi - t£, bj - «* ] U [a, + t£, &t + «*, ^ = (jjja;, 6|]J 

is also locally unbiased and its covariance matrix is, 

V[M'] = V[M] + V . 



□ 



C proof of lemma 7 

Proof The equation ( pl| ) and the equation fl28| ) implies that, for any 
element V of V, there isamx 2m matrix [/ which satisfies 

y- 1 = V _1 (17) = (ReC/L) T Ref/L, 

and 

£/*£/ = Im. (59) 
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Because the map V _1 (*) is continuous and the totality of the m x 2m matrix 
U satisfying (|5^) is compact, the region of V -1 ([7) is compact. Therefore, 
the intersection of V and the set 

{A\A< V }. 

is compact for any real symmetric matrix Vo, for the map 

v- 1 -> V 

is continuous on the intersection of the region of V~ l {U) and the set 

{A\A>V, 1 }, 

both of which are compact. Because Vo is an arbitrary real symmetric 
matrix, we have the lemma. □ 
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